Temperature variation of the pseudogap in underdoped cuprates 
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We consider thermal evolution of the spectral function A{Q,T) in underdoped cuprates at T > Tc- 
We find that in the strong coupling limit, the fermionic Green's function near (0, tt) behaves as 
G^^(Jl) oc \/Tl in the whole frequency range relevant to photoemission experiments. The analysis of 
the pairing problem with this form of G yields a leading edge gap in A{Q,) and a broad maximum at 
larger frequencies. We find that classical fluctuations predominantly add a frequency independent 
term to G~^(f2). This term destroys the leading edge gap at T = r*(k) which has a maximum 
value of about 200-?^ near (0, tt). 



The pseudogap phenomenon in the normal state prop- 
erties of underdoped cuprates is one of the key issues in 
the field of high-Tc superconductivity. NMR, transport, 
optical and photoemission measurements ||l| all demon- 
strate that underdoped cuprates behave in many respects 
like a d— wave superconductor at temperatures which can 
be several times larger than actual Tc. In photoemission 
measurements, the evidence for a pseudogap is based on 
the observation of a leading edge gap (LEG) in the quasi- 
particle spectral function A{Q) near (0, tt) and symmetry 
related points [^|J^. Contrary to situation below Tc, how- 
ever, the LEG is not accompanied by a sharp quasipar- 
ticle peak above it. Instead, the spectral function is flat 
above the LEG and only displays a broad maximum at 
much larger frequencies [1). This form of A{Q) exists be- 
tween Tc, which can be very small in strongly underdoped 
cuprates, and some higher temperature, T*{k) which is 
momentum dependent and has a maximum value of 
about 200K near (0,7r). Similar observations have been 
made recently in tunneling experiments [Q. 

In a recent paper ^ , one of us argued that the LEG 
naturally appears in the in the "magnetic" scenario for 
cuprates due to strong quantum magnetic fluctuations 
which give rise to d— wave precursors above Tc, but si- 
multaneously almost destroy the quasiparticle peak near 
(0,7r). These two effects yield a LEG without a quasi- 
particle peak, in agreement with the data. 

In this paper we show that, as temperature increases, 
magnetic fluctuations loose their dynamics and progres- 
sively acquire a classical, quasistatic behavior. We show 
that this effect acts against d— wave precursors. We fur- 
ther show that the LEG is filled in by excitations at 
temperatures which are comparable to a typical spin- 
fluctuation frequency and also depend on the momentum, 
being largest near (0,7r). 

As in earher studies (5|-|^], our consideration is based 
on the spin-fermion model in which itinerant fermions 
which form a Fermi liquid with a large, Luttinger-type 
Fermi surface interact with their collective spin degrees 
of freedom by 
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Here g is the coupling constant which is assumed to 
increase as the system approaches half-filling, Oi are 
the Pauli matrices, and S'q are collective spin variables 
which are characterized by their full spin susceptibility 
X^5(g,t^) = ^^iXl'y^ + (qC)^ - ^t^Mf - u?!^^). Here ^ 
is the correlation length, A cx and q = Q — q where 
Q is the momentum at which the static susceptibility 
is peaked. At small doping, Q is either equal or very 
close to the antiferromagnetic momentum (tt, tt). Fur- 
ther, X = Q-C^ is the static susceptibility at q = Q and 
Wsf is the typical frequency of overdamped spin fiuctua- 
tions. A fit to NMR experiments yields a ~ 25 states/eV, 
ciJsf ^10 meV and ^ ~ 3 — 4 at the onset temperature for 
a LEG 1^ . The value of A is unknown, but since term 
in X accounts for spin waves, it is reasonable to estimate 
A ~ WisC^^ where Vs is the spin- wave velocity. 

Our goal is to obtain the form of the quasiparticle 
spectral function ^(k, fi) = (l/7r)|Im G(k, ri)| where 
G'-i(k,17) = Go^(k,r2) - S(k,17). This can be most 
straightforwardly done by first evaluating E(k, fim) for 
discrete Matsubara frequencies = 7rr(2m -I- 1), and 
then performing analytical continuation to real frequen- 
cies. To avoid unnecessary complications, we also re- 
strict our consideration to the points in momentum space 
k = khs where ek = Ck-i-Q — Mi ^iid the self-energy has 
the largest value. In cuprates these points (hot spots) are 
located near (0, tt) and symmetry related points. Finally, 
we restrict our consideration to second-order term in g 
- previous studies indicate that the second-order result 
for the fermionic self-energy is dominant even for large 
values of the coupling constant 

To second order in g, S(k, r2m) is given by 
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where w„ = limT and Go is a bare fermionic Green's 
function which near the Fermi surface has a form 
Go(k,a;„) = Zo/(icj,i - ek). Expanding in ek+q to lin- 
ear order in a deviation from a hot spot and performing 
integration over momentum we obtain 
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where = l+|w„|/wsf +u;^/a2_((J7„+^„)/(„^-1))2^ 
and we introduced (7q = ^■kv£,~^ / (iaZ^) where v is the 
Fermi velocity at a hot spot. The fits to photoemission 
and NMR data above T* yield, respectively, v ~ OAeV 
and v^~^/oJs{ ~ 10 independent on temperature. For 
definiteness, we will also choose v^^^/A = 1. We have 
done calculations for larger ratios (2 and 4) and obtain 
little difference in the results. 

The summation over n can be performed explicitly for 
rim ^ v(,~^. For these ftm, the use of the Poisson sum- 
mation formula yields 

where A(f2m, T) is expressed in terms of Fresnel integrals. 
The first term in this formula is a quantum result while 
the second term accounts for thermal effects. At low T < 
LJsi, the use of the asymptotic form of Fresnel integrals 
yields 

A(a„, T) = signf^„, ^J^/ (^) , (3) 

where /(O) = 1 - 0.32(7rr/w,f)2 + .... and f{x > 1) = 
(2/3)(l-0.28(7rr/wsf)2 + ... . The leading term in the ex- 
pansion obviously yields a conventional Fermi-liquid re- 
sult for the self-energy. We see, however, that the typical 
scale above which Fermi-liquid expansion breaks down 
for both external frequency and temperature is given by 
ujs{ which is rather small in cuprates. In the opposite 
limit of "high" T, r2 > uj^i (but stiU, r,r2 < the 
self-energy has the form 

E(khs,f^) = -2 (^^^ (r!^3f)'/'(l + »a(^^)), (4) 

where a{n) = 1 + P^/uj^i/U, and p = 7rT/(2wsf)- In 
this limit, the spin fluctuations are essentially static, the 
thermal part of S(khs, ^) dominates, and the self-energy 
is linear in T. We however studied the crossover between 
quadratic and linear in T behavior of f3 in more detail 
and found that although the deviations from the Fermi- 
liquid behavior begin at rather small T, the asymptotic, 
linear in T form of (3 is recovered only at T ^ ajgf, i.e, 
the crossover region is numerically rather wide. 

We will show below that the onset temperature for 
the LEG falls into the intermediate temperature range 
where quantum and classical contributions to the self- 
energy are of about the same strength. But before, we 
briefly review the earlier results of a purely quantum 
consideration when (3 — which was argued to give rise 



to a LEG without a quasiparticle peak. The key points 
of the quantum consideration are: 

(i) The self-energy as in (Q) survives in some region 
around a hot spot with the width in k-space ~ fl/ujsi^'^- 
Numerical estimates show that for > Wsf, this region 
is rather wide and, e.g., for Q = (7r,7r) it includes (0,7r) 
and symmetry related points. 

(ii) Whcn g becomes few times larger than 170 (which 
is likely to happen when the system approaches opti- 
mal doping from the overdoped side), the self-energy 
overshadows the bare fermionic dispersion such that 
G~^{'k,fl) « S(khs,f^)- As real and imaginary parts of 
S(khs,fi) in (^) are the same for /3 = 0, the dominance 
of the self-energy over bare dispersion leads to a trans- 
form of a spectral weight to higher frequencies, and to a 
near-destruction of the quasiparticle peak near hot spots. 

(iii) The analysis of the pairing susceptibility shows 
that though the d— wave pairing susceptibility turns out 
to be substantially enhanced for g > go, the near- 
destruction of the quasiparticle peak prevents the sys- 
tem to become an actual superconductor down to very 
small T. At the same time, due to the enhancement and 
narrowing of the pairing susceptibility at g > 170, the 
frequency shift which one obtains by solving the Gorkov 
equation for the full Green's function in the pairing chan- 
nel becomes larger than the width of the pairing suscep- 
tibility. In this situation, the solution of the Gorkov's 
equation has almost the same form as in the actual su- 
perconducting state where the pairing susceptibility has 
a (5— functional form, i.e., system develops precursors to 
the d— wave superconducting state without developing a 
true superconductivity. It is essential, however, that as 
G^^(k, fl) cx e*'^/'*\/f2, the frequency shift due to c?— wave 
precursors yields a pole in the full G along an imaginary 
rather than a real frequency axis. This in turn yields the 
spectral function in the form 

^ r2 + fin 

Ap=o{n) « Vn ^^^^ (5) 

where explicit computations yielded ||^ fio = ^oik) ~ 
0.68a;sf^^(cos fca; — cosfcj,)^. This spectral function does 
not possess a sharp quasiparticle pole, but nevertheless 
it shows a "gap-like" behavior. Indeed, A{uj) rapidly 
increases at low frequencies, reaches half maximum at 
~ O.2f2oj then passes through a broad maximum at 
= r^o, and slowly decreases at higher frequencies. This 
behavior, shown in Fig. 0, is in a good agreement with 
the form of the photoemission curve observed in experi- 
ments below T* d. 

We now study how the spectral function changes as 
one takes thermal effects into consideration. We assume 
that the conditions for c?— wave precursors are satisfied 
(i.e., the width of the pairing susceptibility is smaller 
than the frequency shift obtained by solving the Gorkov's 
equation), but the self-energy, and, hence inverse bare 
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FIG. 1. The plot of the spectral function A{Q) given by 
Eq. Solid, dashed-dotted and dotted lines correspond 

to /? = 0, 1/2 and 1, respectively. The spectral function for 
/3 = displays the LEG and the broad maximum at higher 
frequencies. As f3 increases with temperature, the LEG is 
progressively filled in by excitations. 
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FIG. 2. The results for the full self-energy (||) for first 
8 Matsubara frequencies i},n ~ ■nT{2m + 1) and the fits 
to \/Q,m + P^JiOsf /2 dependence. Solid, dashed-dotted and 
dashed lines correspond to T = ciJsf/2, T — ojsf and T — 2uisi, 
respectively. 



Green's function for the pairing channel, has the form of 
Eq. (|^) with P > 0. For the spectral function, we then 
obtain instead of (||) 

Here it is convenient to rewrite a{^l) as a(0) = 1 + 
f3^/no/n where f3 = l3/{0.82^\cosk^ - cosky\). The plots 
of A{fl) for various values of P are presented in Fig. |^. We 
see that the LEG is progressively filled in as f3 increases 
from 0.5 to 1. This happens chiefly because thermal fluc- 
tuations eliminate the a/TJ dependence at small frequen- 
cies. Glearly, the temperature where the LEG disappears 
depends on k being the largest near the points where the 
d— wave gap has the largest value. The high-frequency 
part of A{Q), which possesses a broad maximum, is a 
bit more robust and disappears at somewhat larger (3. 
In other words, as temperature increases, thermal exci- 
tations first fill the LEG, and then, at somewhat higher 
T, destroy the broad maximum at Qq. This behavior 
seems to be consistent with both photoemission and 
tunneling experiments Q . Notice that this mechanism is 
different from a conventional thermal smearing of the gap 
which in our approach would would correspond to a ther- 
mal broadening of the pairing susceptibility, x"^{^)- We 
now need to find the actual temperature dependence of 
(3 and obtain the onset temperature for the LEG which 
can be directly compared with the data. For this, we 
compute the self-energy numerically by performing the 
summation over n in without simplifying the form of 
the integrand, i.e., without assuming that Vim <C w^"^. 
We performed computations for T = Wgf/SjT = Wgf and 
T = 2usf and in all three cases found that the self-energy 



for external fermionic frequencies f2„j — 7rT(2m + 1) and 
< TO < 7 can be well approximated by the same form 
as in (I), i.e. S(khs,rJ™) cx (V^U + /3Jujsf/2. The 
accuracy of the fit is demonstrated in Fig. g. We empha- 
size that this procedure shows that the exact self-energy 
follows the \/f2 dependence in a much larger frequency 
range than one may expect from the approximate analyt- 
ical treatment. The overall factor in the fit decreases a 
bit with increasing temperature. This last renormaliza- 
tion, however, can be absorbed into the renormalization 
of g and is irrelevant as the positions of both the LEG 
and the broad maximum in AiTl) are governed by fioj 
which does not depend on the coupling strength as long 
as g> go. 

We now present the numbers. We obtained f3 = 
0.24(0.38), 1.69(2.57) and 6.82(10.63) for ^-^/A = 1(2) 
and T/uJsf = 0.5,1 and 2, respectively. We see that (3 
rapidly increases with increasing temperature. The in- 
crease roughly follows T'^ dependence with v somewhat 
larger than 2, i.e., it is much stronger than the linear in 
T increase which one would expect for a purely classical 
self-energy. On the other hand, already for T — uJsf/2, 
the value of /? is about 40% smaller than in the Fermi- 
liquid theory. This collaborates our earlier observation 
in (^) that the Fermi-liquid behavior is altered already 
at low T, but there is a wide intermediate temperature 
range where the thermal part of the self-energy interpo- 
lates between quadratic and linear in T behavior. 

Using f3 = 0.5 as a criteria for the onset of the 
LEG, we find that the LEG near (0, tt) is filled in at 
T - T* = 0.6wstC^/''- A fit to NMR shows that 
T* - 150 - 200K. Away from (0,7r), the LEG is 
filled in at smaller temperatures which roughly follows 
T*(k) = T*{\cosk^ - cos/cy 1/2)1/''. Both of these re- 
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suits are in quantitative agreement agreement with the 
data H,^. In addition, as one moves towards the zone 
diagonal, i.e., goes away from the hot spots, the over- 
all strength of the term in the self-energy decreases. 
This pushes the system back towards a weak coupling 
limit and causes an additional decrease in T*(k). Simul- 
taneously, one should gradually recover some features of a 
conventional quasiparticle peak in A(Q). This last effect 
has also been observed in recent experiments . Finally, 
a fit to the NMR data shows that the ratio T/cUsf flat- 
tens at ~ 1.4—1.5 slightly above T* and remains constant 
upto much higher T. Then, P also flattens at a value ^ 1. 
This /3 still yields a broad maximum in A{n) above T*, 
in agreement with the data. 

We also emphasize that in our approach, the LEG 
gradually appears at T ^ T*(k) without a break up of 
the Fermi surface. In other words, the present theory 
predicts that above Tc, there still exists a quasiparticle 
peak at a vanishing frequency at k = kp. The residue 
of the peak, however, decreases as one moves away from 
the zone diagonal. Near (7r,0) it scales as (ffo/ff)^ which 
seems too small to be detected by current photoemission 
experiments. Notice in this regard that photoemission 
data show that the quasiparticle peak is clearly present 
near the zone diagonal, and gradually disappears as one 
approaches (0,7r). The tunneling experiments which 
measure the momentum averaged spectral density, also 
seem to observe a coherent peak above Tc in the unoccu- 
pied part of the density of states. 

A final note. In the present consideration, the LEG 
and the broad maximum in A(uj) at a finite frequency 
are both d— wave precursors in a Fermi liquid with strong 
magnetic fluctuations. A number of other scenarios for a 
pseudogap have also been proposed [pPj9|,p^ ■ In partic- 
ular, several groups argued that the broad maxi- 
mum in the spectral function is a precursor of an anti- 
ferromagnetic SDW state. In the present calculations, 
SDW precursors do not appear because the fermionic 
self-energy is an increasing function of frequency upto 
very large frequencies. The SDW precursors, on the con- 
trary, can appear only if the self-energy decreases with 
frequency j^. We found, however, that the form of the 
spectral function at > v^~^ is highly sensitive to the 
high-frequency form of the spin susceptibility. Namely, 
we did calculations using xCq? ^) in which the term is 
substituted by a hard cutoff at w ~ A and found that the 
classical behavior sets up at much smaller temperatures. 
In this last case, for T > T* , the dominant contribution 
to comes from the n = term, i.e., spin fluctua- 
tions are essentially static. These fluctuations yield the 
fermionic self-energy which decreases with frequency and 
gives rise to SDW precursors at > w^^^ provided that 
i.9/9of{T/v£,-'^) > 1/4 which is satisfied above T* §. 
Quantum SDW precursors can also appear at high fre- 
quencies provided that feedback effects substantially re- 
duce the spin damping |^,^,^. In view of this nonuni- 



versality, it is possible that the broad maximum in A{uj) 
at high frequencies is actually a mixture of d-wave and 
SDW precursors. 

To summarize, in this paper we discussed the ther- 
mal variation of the fermionic spectral function A{Q) in 
the spin-fermion model for underdoped cuprates. We 
found that in the strong coupling limit, the fermionic 
Green's function near (0, tt) behaves as G~^{il.) (x (\/f2 + 
(3y^uJsf in the whole frequency range relevant to 
photoemission experiments. The frequency-independent 
term is due to thermal spin fluctuations. At low temper- 
atures, the Vfi term dominates and yields the LEG in 
A{fl) and a broad maximum at larger frequencies. As T 
increases, spin fluctuations progressively acquire a clas- 
sical, quasistatic behavior. This effect yields a rapid in- 
crease of (3 with temperature which in turn gives rise to 
a destruction of the LEG. We found that the LEG disap- 
pears at T = T* - 0.6cJsfC^/'' which is 150 - 200K near 
(0, tt) and decreases as one moves towards zone diagonal. 
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